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I. INTRODUCTION 

Real atoms are typically complex, having ground and 
excited states with spin structure. The molecules formed 
from the atoms typically have a rich spectrum of near- 
threshold bound and quasi-bound molecular states when 
the molecular spin, rotational, and vibrational structure 
is taken into account. When an ultracold gas of atoms 
is produced, the atoms are prepared in specific quantum 
states, and collisions between the atoms occur with an ex- 
tremely precisely defined energy close to the E = colli- 
sion threshold of the interacting atoms, where E denotes 
energy. The collision then makes the near-threshold spec- 
trum of the molecular complex of the two atoms accessi- 
ble to electromagnetic probing. An external magnetic or 
electromagnetic field can be precisely tuned to couple the 
colliding atoms to a specific molecular state, which can be 
viewed as a scattering resonance. This permits both ex- 
traordinary spectral accuracy in probing near-threshold 
level positions (order of E/h = 10 kHz accuracy for 1 /iK 
atoms) and precise resonant control of the collisions that 
determine both static and dynamical macroscopic prop- 
erties of quantum gases. Consequently, understanding 
the near-threshold bound and scattering states is essen- 
tial for understanding the collisions and interactions of 
ultracold atoms. This is also true for interactions of ul- 
tracold molecules. 

This Chapter concentrates on understanding molecules 
that can be made by combining two cold atoms using ei- 
ther magnetically tunable Feshbach resonance states \M 
or optically tunable photoassociation resonance states [2j • 
Such resonances provide a mechanism for the formation 
of ultracold molecules from already cold atoms. In ad- 
dition, magnetically tunable resonances have been used 
very successfully to control the properties of ultracold 
quantum gases. This Chapter treats both magnetically 
and optically tunable molecular resonances with the same 
scattering theory framework. The viewpoint from quan- 
tum defect theory is emphasized of conceptually sepa- 
rating the interaction of the atoms into short range and 
long range regions. These regions are characterized by 
very different energy and length scales. Much insight 
about near-threshold collisions and bound states, as well 
as practical tools for their study, can be gained b y ta king 
advantage of this separation [|, i, H, H, 0, S, i,M EH- 
While molecular physics is typically concerned with 
strong short range interactions associated with " ordinary 
molecules," ultracold physics is concerned with scatter- 



ing states and very weakly bound molecular states in the 
threshold domain near E = 0. The long range poten- 
tial, which has a lead term that varies as l/R n , plays an 
important role in connecting these two regimes. 

We briefly summarize here the theory of cold collisions, 
which is described in detail in Chapter XXX. The scat- 
tering wavefunction is expanded in states of relative an- 
gular momentum of the two atoms characterized by par- 
tial wave quantum number £ = 0, 1, 2 . . .. Generally, the 
atoms can be initially prepared in one of several quan- 
tum states, and the scattering "channels" can be speci- 
fied by a collective set of quantum numbers a represent- 
ing the state of each atom and the partial wave. Upon 
solving the Schrodinger equation for the system, the ef- 
fect of all short-range interactions during a collision with 
E > is summarized in the scattering wavefunction for 
R — > oo by a unitary S'-matrix. Only the lowest few 
partial waves can contribute to cold collisions, and in 
the limit E — > 0, only s-wave channels with £ = have 
non-negligible collision cross sections. Using the complex 
scattering length a — ib to represent the s-wave S'-matrix 
element S aa = exp [— 2ik(a — ib)] in the limit E — > 0, the 
contribution to the elastic scattering cross section from 
s-wave collisions in channel a is 

o-el = km 9^n\l- S aa \ 2 = Agir (a 2 + b 2 ) , (1) 



where Tik = ^/2/iE is the relative collision momentum in 
the center of mass frame for the atom pair with reduced 
mass \i. The rate coefficient A"i oss = u\ a&& v for E — > s- 
wave inelastic collisions that remove atoms from channel 
a is 

K losB = lim g^ (1 - \S aa \ 2 ) = 2. 9 -6 (2) 

-E^O fik [l 

where v = Hk/n is the relative collision velocity. The 
symmetry factor g = 1 when the atoms are bosons or 
fcrmions that are not in identical states, g = 2 or g = 1 
respectively for two bosons in identical states in a normal 
thermal gas or a Bosc-Einstcin condensate, and g = for 
two fermions in identical states. If there are no exocrgic 
inelastic channels present, then 6 = and only elastic 
collisions are possible. 

The Schroodinger equation also determines the bound 
states with discrete energies Ei < 0. While the con- 
ventional picture of molecules counts the bound states 
by vibrational quantum number v = 0, 1 . . . from the 
lowest energy ground state up, it is more helpful for the 
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present discussion to count the near-threshold levels from 
the E = dissociation limit down by quantum numbers 
i = — 1, —2 . . .. In the special case where a — > +oo, the 
energy of the last bound s-wave state of the system with 
i = — 1 depends only on a and \x and takes on the follow- 
ing "universal" form: 



2fia 2 



(3) 



Section |TT] describes the bound and scattering prop- 
erties of a single potential with a van der Waals long 
range form. Section IIIII extends the treatment to multi- 
ple states and scattering resonances. Sections [IV] and [V] 
respectively discuss the properties of magnetically and 
optically tunable molecular resonance states. 



II. PROPERTIES FOR A SINGLE POTENTIAL 

In this section let us ignore any complex internal 
atomic structure and first consider two atoms A and B 
that interact by a single adiabatic Born-Oppenheimer 
interaction potential V(R), illustrated schematically in 
Fig. [TJ The wavefunction for the system is \a)\ipt) / R, 
where \a) represents the electronic and rotational degrees 
of freedom, and the wavefunction for relative motion is 
found from the radial Schrodingcr equation 
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Solving Eq. |4] gives the spectrum of bound molecular 
states tpu with energy En = — h 2 k 2 ( /(2^) < and 
the scattering states ipe(E) with collision kinetic energy 
E = h 2 k 2 /(2[i) > 0, where ku and k have units of 
(length) -1 . As R — > oo, the bound states decay as e~ kaR 
and the scattering states approach 



ipt(E) -> csin(fc J R-7r^/2 + ?7£)/fc 1/2 . 
Bound states arc normalized to unity, \(il>u\ipj£' 



(5) 



We choose the normalization constant c = 



2h/Tl 2 it so that scattering states are normalized per 
unit energy, (ip e (E)\ip e (E')) = 5(E - E')5 U '- Thus, the 
energy density of states is included in the wavefunction 
when taking matrix elements involving scattering states. 

The long range potential between the two atoms varies 
as —C n /R n . We are especially interested in the case of 
n = 6 for the van der Waals interaction between two 
neutral atoms. This is the lead term in the long-range 
expansion of the potential in inverse powers of R that 
applies to many atoms that are used in ultracold exper- 
iments. This potential has a characteristic length scale 

of i? v dw = y 2/j,Ce/h 2 /2 that depends only on the values 

of fi and Cq Q. Values of Cq are tabulated by Dere- 
vianko [l2| for alkali-metal species and by Porsev and 
Derevianko [l3[ for alkaline-earth species. We prefer to 



TABLE I: Characteristic van der Waals scales a and E for 
several atomic species. (1 amu =1/12 mass of a 12 C atom, 1 
au= 1 Eh,a% where Eh is a hartree and 1 ag— 0.0529177 nm) 





lllclSS 


c 6 


a 


E /h 






(amu) 


(au) 


(oo) 


(MHz) 


(mK) 


6 Li 


6.015122 


1393 


29.88 


671.9 


32.25 


23 Na 


22.989768 


1556 


42.95 


85.10 


4.084 


40 K 


39.963999 


3897 62.04 


23.46 


1.126 


87 Rb 


86.909187 


4691 
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6.668 


0.3200 


88 Sr 
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0.3827 


133 Cs 


132.905429 
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96.51 


2.916 


0.1399 


174 Yb 


173.938862 


1932 


75.20 


3.670 


0.1761 
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FIG. 1: Schematic figure of the potential energy curve V(R) 
as a function of the separation R between two atoms A and 
B. The horizontal lines labeled AB indicate a spectrum of 
molecular bound states leading up to the molecular dissocia- 
tion limit at E — 0, indicated by the dashed line. The long 
range potential varies as ~C n /R n ■ 



use a closely related van der Waals length introduced by 
Gribakin and Flambaum \\M 



= 47r/T(l/4) 2 i? vdw = 0.955978 ...R 



vdw ■ 



(0) 



where T(x) is the Gamma function. This length defines a 
corresponding energy scale E = h 2 /(2/ia 2 ). The param- 
eters a and E occur frequently in formulas based on the 
van der Waals potential. The wavefunction approaches 
its asymptotic form when R 3> a and is strongly influ- 
enced by the potential when R < a. Tabic [T] gives the 
values of a and E for several species used in ultracold 
experiments. 

Samples of cold atoms can be prepared with kinetic 
temperatures on the order of nK to mK. The energy asso- 
ciated with temperature T is ksT where ks is the Boltz- 
mann constant. For example, at T = 1 /iK, ksT = 0.86 
neV and ksT/h = 21 kHz. This ultracold energy scale 
is 9 to 10 orders of magnitude smaller than the en- 
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FIG. 2: Radial wavefunction ip (R) for I = at E/k b = 1 
/xK for the pairs 174 Yb- 174 Yb (solid), 171 Yb- 171 Yb (dashed), 
and 170 Yb- 173 Yb (dotted), which have respective scattering 
lengths of 105 ao, -3 ao, and —81 ao [l5T |. The inset shows an 
expanded view of the wavefunction on a smaller length scale 
on the order of a, the characteristic length of the van der 
Waals potential. The 174 Yb- 174 Yb case shows the oscillations 
that develop when R < a. 



FIG. 3: wavefunctions for the last s-wave i = — 1 bound state 
(solid line) with E-i t o/h = —10.6 MHz and for the s-wave 
scattering state (dashed line) for E/h = 0.02 MHz (E/k B = 1 
/iK) for two 174 Yb atoms. Both wavefunctions are given a 
common JWKB normalization at small R <a and are nearly 
indistinguishable for R < a. The potential supports N — 72 
bound states, and the wavefunction for this i = — 1 and v = 71 
level has N — 1 = 71 nodes. 



ergy scale of 1 to 10 eV associated with ground or ex- 
cited state interaction energies when a molecule is formed 
at small interatomic separation i?bond on the order of 
a chemical bond length. In a cold collision, the ini- 
tially separated atoms have very low collision energy 
E = h 2 k 2 /(2[i) ?s and very long de Broglie wave- 
length 27r/ fc. The atoms come together from large dis- 
tance R and are accelerated by the interatomic potential 
V(R), so that when they reach distances on the order of 
-Rbond they have very high kinetic energy on the order of 
|V(i2 b ond)|- The local de Broglie wavelength 2it/k(R, E) 
in the short range classical part of the potential, where 
k(R,E) = y/2fi(E- V{R))/h, is orders of magnitude 
smaller than the separated atom de Broglie wavelength 
and is nearly independent of the value of E, which is close 
to 0. 

This separation of scales is illustrated in Fig. [21 which 
shows examples of s-wavefunctions at a collision energy 
E/ks = 1 AtK, where dividing E by fcs allows us to ex- 
press energy in temperature units. This example uses 
three isotopic combinations of pairs of Yb atoms, which 
has a spinless electronic configuration and a sin- 
gle ground state electronic Born-Oppcnhcimer potential 
V(R). The species Yb makes a good example case to il- 
lustrate the principles in this section, since it has 7 stable 
isotopes and 28 different atom pairs of different isotopic 
composition for which the threshold properties have been 
worked out All combinations have the same V(R) 

but different reduced masses. This mass-scaling approx- 
imation, which ignores very small mass-dependent cor- 
rections to the potential, is normally quite good except 



for very light species such as Li. Fig. [2] shows that the 
three examples have similar phase-shifted sine waves with 
a common long de Broglie wavelength of 2n/k = 6300 
ao. For small R where kR <C 1 the sine function van- 
ishes as csinfc(i? — a)/v« — ► c^/k(R — a). The actual 
wavefunction oscillates rapidly at small R due to the in- 
fluence of the potential. Since the asymptotic form for 
kR <C 1 varies as fc 1 / 2 as k — > 0, the short range oscil- 
lating part also has an amplitude proportional to fc 1 / 2 
in order to connect smoothly to the asymptotic form as 
k — > 0. This property ensures that the threshold ma- 
trix elements that characterize Feshbach resonances and 
s-wave inelastic scattering are proportional to k 1 / 2 . 

Figure[3]illustrates more clearly the nature of threshold 
short range scattering and bound state wavefunctions. 
When given an appropriate short range normalization, 
near-threshold scattering and bound state wavefunctions 
have a common amplitude and phase in the region of R 
small compared to the range a of the long range potential. 
While this can be put on a rigorous quantitative ground 
within the framework of quantum defect theory Q, it 
is easy to show using the familiar JWKB approxima- 
tion @, S 0| ■ We can always write the wavefunction in 
phase-amplitude form ipi(R,E) = ae(R, E) sin f3e(R, E) 
and transform the Schrodinger equation ^ into a set 
of equations for a>e and fy. The asymptotic ipi(R,E) in 
Eq. ([5]) clearly corresponds to this form with ag — > c/k 1 ^ 2 
as R — > oo. Another familiar form is the JWKB semi- 
classical wavefunction ip^ WKB (R, E), for which 

aj WKB (R,E) = c/UR.E) 1 ' 2 (7) 
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(R,E) 



k e (R',E)dR' + -. (8) 

Rt 4 



where Rt is the inner classical turning point of the po- 
tential. 

When the collision energy E is sufficiently large, so 
there are no threshold effects, the JWKB approximation 
is a excellent approximation at all i?, and the form of 
aj WKB (R 1 E) in Eq. (J7J) applies at all R, transforming 
into the correct quantum limit as R — > oo. On the other 
hand, the JWKB approximation fails for s-waves with 
very low collision energy. This failure occurs in a region 
of R near a and for collision energies E on the order of E 
or less. The consequence is that the JWKB wavefunction, 
with the normalization in Eq. ((7|) , is related to the actual 
wavefunction, with the asymptotic form in Eq. by a 
multiplicative factor Cg{E) : so that as E — > 



MR,E) = C t {E)- l ^i WKB {R,Q) . 



(9) 



As k — > for a van der Waals potential varying as l/R 6 , 
the s-wave threshold form is Co(E)~ 2 = ka[l + (r — l) 2 ], 
where r = a/a is the dimcnsionless scattering length in 
units of a @ . Equation © gives an excellent approxima- 
tion for the threshold ipo(R,E) for R < a and k < 1/a. 
At high energy, when E 3> E, Co(E)^ 1 approaches unity 
and the JWKB approximation for ipo(R, E) applies at all 
R. 

The unit normalized bound state wavefunction ij)n(R) 
can be converted to an " energy normalized" form by mul- 
tiplying by Idi/dEu] 1 / 2 , where —di/dEu > is the en- 
ergy density of states. Away from threshold, this is just 
the inverse of the mean spacing between levels, whereas 
for s-wave levels near threshold for a van der Waals po- 
tential , di/dE. ia — > r/(2nE)~ 1 as fc_i,o = 1/a — > [1]. 
The relation of ipu to the energy- normalized JWKB form 
in the classically allowed region of the potential is 



ipu(R,Etf) 



Oi 



dE lt 



-1/2 



(10) 



Figure [3] plots C (E)ip (R, E) « ip$ WKB (R,0) for 
the scattering state and \di/ dE^^Tp^R, E m ) tn 
tpi WKB (R, 0) for the i = -1 bound state. Thus the near- 
threshold bound and scattering wavefunctions, when 
given a common short range normalization, are nearly 
identical and are well approximated by ipQ WKB (R, 0) in 
the region R < a. For R > a the wavefunctions begin to 
take on their asymptotic form as R — > oo. The shape of 
the wavefunction at very small R on the order of -Rbond is 
usually independent of E for ranges of E /ks on the order 
of many K. The short range shape is even independent 
of £ for small t, since the rotational energy is very small 
compared to typical values of T^(-Rbond)- However, the 
amplitudes of the wavefunctions depend strongly on the 
whole potential, which determines a, and are analytically 
related to the form of the long range potential. 

The separation of scales for R > a and R < a is a key 
feature of ultracold physics that enables much physical 



insight as well as practical approximations to be devel- 
oped about molecular bound and quasibound states and 
collisions. Given that Cq, [a, and the s-wave scattering 
length a are known, the Schrodinger equation ^ can be 
integrated inward using the form of Eq. ([5]) as k — > as 
a boundary condition, thus giving the wavefunction and 
nodal pattern for R < a as E — ► 0. Assume that it is pos- 
sible to pick some R = R m such that i?bond <C R m *C a 
and V(R m ) is well-represented by its van der Waals form. 
Then the log of the derivative of the wavefunction at R m , 
which also can be calculated, provides an inner boundary 
condition, independent of E over a wide range of E, for 
matching the wavefunction at E propagated from large 
R. All that is needed to do this is to know a, /i and 
the long range potential. Thus, it is readily seen that all 
of the near threshold bound and scattering states, even 
those for £ > 0, can be calculated to a very good approx- 
imation for R > R m once Cq, /i, and a are known. 

Figure [4] shows the spectrum of bound states E^, in 
units of E, for £ up to 5 for two cases of scattering length, 
based on the van der Waals quantum defect theory of 
Gao 0, [HJ]. Panel (a) shows the case of a = ±oo, where 
there is a bound state at E = 0. The locations of the 
bound states for a = ±oo define the boundaries of the 
"bins" in which, for any a, there will be one and only one 
s-wave bound state, for example, —36. IS < S_i q < 
and -249E < E_ 2 ,o < -36.125. The panel also shows 
the rotational progressions for each level as £ increases. 
The a = ±oo van der Waals case also follows a "rule of 4", 
where partial waves £ — 4, 8, . . . also have a bound state 
at E = 0. Panel (b) shows how the spectrum changes 
when a = a, for which the there is a d-wave level at 
E = 0. Similar spectra can be calculated for any a. 

Gribakin and Flambaum [3| showed that the near- 
threshold s-wave bound state for a van der Waals po- 
tential in the limit a > a is modified from the universal 
form in Eq. (J3j) as 



E. 



2/i(a — a)' 



(11) 



This approaches the universality limit when a > d, in 
which case the s-wave wavefunction takes on the univer- 
sal form i/j (R,E) = y/2/ae~ R / a . Such an exotic bound 
state, known as a "halo molecule," exists primarily in 
the nonclassical domain beyond the outer classical turn- 
ing point of the long-range potential with an expecta- 
tion value of R of a/2, which grows without bound as 
a -> +oo [J]. 

Bound state and scattering properties are closely re- 
lated. It is instructive to imagine that there is some 
control parameter A that can be varied to make the scat- 
tering length vary over its whole range from +oo and 
— oo, changing the corresponding bound state spectrum. 
One way to do this would be to vary the reduced mass. 
Of course, this is not physically possible. However, there 
are elements with many isotopes, so that a wide range of 
discrete reduced masses are possible. An excellent physi- 
cal system to illustrate this is the Ytterbium atom, used 
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FIG. 4: Dimensionless bound state energies Eu/E for partial 
waves I = 0...5 (s,p, d, /, g, h). Panel (a) is for the case 
a = ±oo and Panel (b) is for a — a. 
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FIG. 5: The upper panel shows s-wave scattering length and 
the lower panel shows bound state binding energies —Eio(X) 
for Yb2 molecular dimers versus the control parameter A = 
2/i. The vertical dashed lines show the points of singularity 
of a(A). The horizontal dashed lines show the boundaries of 
the bins in which the i = — 1 and i = —2 levels must lie. 



in the examples of Figs. [2] and [3] The stable isotopes with 
masses 168, 170, 172, 174, and 176 are all spinless bosons 
and the 171 and 173 isotopes are fermions with spin 1/2 
and 5/2 respectively. Yb atoms can be cooled into the 
pK domain and all isotopes, including the fermionic ones 
in different spin states, have s-wave interactions.. The lo- 
cations of several £ = and 2 threshold bound states of 
different isotopic combinations of Yb atoms in Yb2 dimcr 
molecules have been measured, and the long range poten- 
tial parameters and scattering lengths determined [TBI ]. 

Figure [5] shows the s-wave scattering length and bound 
state binding energies versus the continuous control pa- 
rameter A = 2/i. Physically, there are 28 discrete values 
between A =168 and 176. The scattering length has a sin- 
gularity, and a new bound state occurs with increasing A, 
at A = 167.3, 172.0, and 177.0. The range between 167.3 
and 172 corresponds to exactly N — 71 bound states in 
the model potential used. Near A = 167.3 the last s- 
wave bound state energy E—± o ^ as — h 2 /(2pa 2 ) as 

a — ► +oo. The binding energy | _E7 x ,o | gets larger as A 

increases and a decreases, so that for a van der Waals 
potential E—\ o approaches the lower edge of its "bin" 
at — 36.1.E as a — > — oo. As A increases beyond 172.0, 
the i — — 1 level becomes the i = — 2 level as a new 
% = — 1 "last" bound state appears in the spectrum. 

The variation of scattering length with 2/i is given by 
a remarkably simple formula. While semiclassical theory 
breaks down at threshold, Gribakin and Flambaum (I4| 
showed that the correct quantum mechanical relation be- 
tween a and the potential is 



a 1 - tan ( $ - — 



(12) 



where 



$ = 



yj -2nV(R)/h 2 = (5 J WKB (oo, 0) - tt/4 . (13) 



The number of bound states in the potential is N — 
[$/7r — 5/8] + 1, where [. . .] means the integer part of 
the expression. These expressions work remarkably well 
in practice. Although the results in Fig. [5] are obtained 
by solving the Schrodinger equation for a realistic po- 
tential, virtually identical results are obtained for a from 
Eq. Q12p. In fact, a and Eio are nearly the same on the 
scale of Fig.[5]if the simple hard-core van der Waals model 
of 14[ is used for the potential, namely V(R) = —Cq/R 6 
if R > Rq and V(R) = +oo if R < R , where the cut- 
off Rq is chosen to fit a or -E_i.o data from two different 
isotopes. With the mass scaling oc ^[Ji in Eq. (Tl3|) . know- 
ing Cq and -E-1,0 for two isotopic pairs determines a and 
£7-1,0 for all isotopic pairs. The approximation is fairly 
good even for levels with larger |z| or I > 0, although it 
will become worse as \i\ or £ increase. 

In summary, it is very useful to take advantage of the 
enormous difference in energy and length scales associ- 
ated with the cold separated atoms and deeply bound 
molecular potentials. This allows us to introduce a gen- 
eralized "quantum defect" approach for understanding 
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threshold physics @, i, 0, 0, EH ■ Threshold bound state 
and scattering properties are determined mainly by the 
long range potential, once the overall effect of the whole 
potential is known through the s-wave scattering length. 
A similar analysis can be developed for other long range 
potential forms, for example, 1/i? 4 ion-induced dipole or 
1/i? 3 dipole-dipole interactions. 



III. INTERACTIONS FOR MULTIPLE 
POTENTIALS 

Generally the cold atoms used in experiments have ad- 
ditional angular momenta (electron orbital and/or elec- 
tron spin and/or nuclear spin), so that more than one 
scattering channel a can be involved in a collision. Each 
channel has a separated atom channel energy E a . Fig. [I] 
could be modified to illustrate such channels by adding 
additional potentials and their corresponding spectra dis- 
sociating to the E a limits. If Etot is the total energy of 
the colliding system, the designation open or closed is 
used for channels with E to t > E a or E to t < E a respec- 
tively. Inelastic collisions from entrance channel a are 
possible to open exit channels (3 when E a > Ep , whereas 
closed channels (3 can support quasibound states as scat- 
tering resonances when E a < E tot < Ep. The ability 
to tune resonance states to control scattering properties 
or to convert them into true molecular bound states is 
an important aspect of ultracold physics that has been 
exploited in a wide variety of experiments with bosonic 
or fermionic atoms [l[. 

Let us first examine the basic magnitude of the s-wave 
inelastic collision rates that are possible when open chan- 
nels are present. The rate constant is determined by the 
magnitude of b in Eq. for which a typical order of 
magnitude is b « a for an allowed transition, that is, 
one with a relatively large short-range interactions in the 
system Hamiltonian. The rate constant can be written 

K loss = 0.84 x 10- 10 5 4^T cm Vs > ( 14 ) 
/z|amu| 

where b is expressed in atomic units (1 au = 0.0529177 
nm) and fj, in atomic mass units (/i = 12 for 12 C). Al- 
lowed processes will typically have the order of magni- 
tude of 10~ 10 cm 3 /s for K\ oss . The s-wave i^i OS s can be 
even larger, with an upper bound of b u = l/(4fc) being 
imposed by the unitarity property of the S'-matrix, i. e., 
< 1 — |5 QCl | 2 < 1. Since the lifetime relative to collision 
loss is t — \/{K\ oss n), where n is the density of the col- 
lision partner, allowed processes result in fast loss with 
t < 1 ms at typical quantum degenerate gas densities. 
This applies to atom-molecule and molecule-molecule col- 
lisions as well as atom-atom collisions. Such losses need 
to be avoided by working with atomic or molecular states 
that do not experience fast loss collisions, such as the 
lowest energy ground state level, which does not have 
exocrgic 2-body exit channels. Alternatively, placing the 



species in a lattice cell that confines a single atom or 
molecule can offer protection against collisional loss. 

An alternative formulation of the collision loss rate is 
possible by rewriting Eq. @ , not taking the E — > limit 
but introducing a thermal average over a Maxwcllian dis- 
tribution of collision energies E, 

^=5^-^E( 1 -| 5 -| 2 )t' (15) 

a 

where Qt is the translational partition function, 1/Qt = 
(27r//fc_BT//i 2 ) 3 / 2 = A?n where At is the molecular ther- 
mal de Broglie wavelength. The (• • -)t expression implies 
a thermal average over the velocity distribution. The sum 
represents a dynamical factor fjj that varies as T 1 / 2 as 
T — > and has an upper bound of unity for s-waves and 
~ ^max if W partial waves contribute at the unitar- 
ity limit. Although Eq. ([15]) reduces to Eq. HU) in the 
T — > s-wavc limit, it lets us see that the collision rate 
is given by an expression having the form 

t- 1 = K loss n = g (nA 3 T )^f D . (16) 

This form embodies some general principles for any colli- 
sions of atoms and molecules. The dimcnsionless nA^ 
factor shows that the collision rate is proportional to 
phase space density of the collision partner (scale by mass 
ratios to convert to an atomic phase space density) . The 
ksT/h factor sets an intrinsic rate scale (dimension of 
inverse time) associated with T. The dimcnsionless fac- 
tor fjj embodies all of the detailed collision dynamics. 
Even using fast time-dependent manipulations to control 
fa does not change the fundamental thermodynamic lim- 
its imposed by the phase space density and ksT/h fac- 
tors. Given Eqs. fT4|) and (| L 5[) and plausible assumptions 
about bor fn, it is possible to estimate the time scales for 
a wide variety of atomic and molecular collision processes 
under various kinds of conditions. 

Now we will examine the important case of tunable 
resonant scattering when a closed channel is present. As- 
sume that open entrance channel a. with E a chosen as 
E a = 0, is coupled through terms in the system Hamil- 
tonian to a closed channel (3 with < E < Ep. Then 
a molecular bound state in channel (3 becomes a quasi- 
bound state that acts as a scattering resonance in chan- 
nel a. Using Fano's form of resonant scattering the- 
ory [l6[, let us assume a "bare" or uncoupled approx- 
imate bound state \C) — ip c (R)\c) with energy E c in 
the closed channel (3 = c and a '"bare" or background 
scattering state \E) = ipb g {R, E)\bg) at energy E in the 
entrance channel a = bg. The scattering phase shift 
rj(E) = rjbg(E) + rj res (E) of the coupled system picks up 
a resonant part due to the Hamiltonian coupling W(R) 
between the "bare" channels. Here r]b g is the phase shift 
due to the uncoupled single background channel, as de- 
scribed in the last Section, and 

n r es(E) = - tan" 1 ( , , (17) 

' v ' \E- E C -SE(E) J 1 v ' 
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has the standard Brcit-Wigncr resonance scattering form. 
The two key features of the resonance are its width 



T(E)=2n\(C\W(R)\E)\\ 



and its shift 



r 

8E{E) =V 

J —i 



A|2 



\(C\W(R)\E ')\ 
E-E' 



■dE' . 



(18) 



(19) 



The primary difference between an "ordinary" reso- 
nance and a threshold one as E — * is that for the for- 
mer we normally make the assumption that T(E) and 
SE(E) are evaluated at E = E c and are independent 
of E across the resonance. By contrast, the explicit en- 
ergy dependence of F(E) and SE(E) are key features of 



threshold resonances 
the E — > limit for s-waves 



1 



111 . 1 1 71 , Il8l | . In the special case of 



-T(E) -» (ka bg )r 
E C + SE(E) -> E , 



(20) 
(21) 



where To and Eq are ^-independent constants. Note that 
T(E) is positive definite, so that r has the same sign as 
ab g - Assuming an entrance channel without inelastic loss, 
so that 7]b g (E) — > —ka.bg, and for the sake of generality, 
adding a decay rate 7 c /7i for the decay of the bound state 
|C) by irreversible loss processes, gives in the limit of 
E-> 0, 



= a — ib = 



a bg 



abg^o 



E Q - »( 7c /2) 



(22) 



This formalism accounts for both kinds of tunable res- 
onances that are used for making cold molecules from 
cold atoms, namely, magnetically or optically tuned res- 
onances. We now give our attention to each of these in 
turn. 



IV. MAGNETICALLY TUNABLE 
RESONANCES 

Cold alkali metal atoms have a variety of magnetically 
tunable resonances that have been exploited in a num- 
ber of experiments to control the properties of ultracold 
quantum gases or to make cold molecules. For the most 
part, experiments have succeeded with species that ei- 
ther do not have inelastic loss channels, or if they do, the 
loss rates are very small. Thus, for practical purposes, 
we can set the resonance decay rate j c = in examining 
a wide class of magnetically tunable resonances. While 
general coupled channel methods can be set up to solve 
the multichannel Schrodinger equation [l|, we will use 
simpler models to explain the basic features of tunable 
Feshbach resonance states. 

Many resonances occur for alkali metal species in their 
2 S electronic ground state because of their complex hy- 
pcrfinc and Zeeman substructure with energy splittings 



very large compared to fcsT. Thus, closed spin channels 
that have bound states near E a of an entrance channel a 
can serve as tunable scattering resonances for threshold 
collisions in that channel. The key to magnetic tuning 
of a resonance is that the resonance state |C) has a dif- 
ferent magnetic moment fi c than the moment /i a toms of 
the pair of separated atoms in the entrance channel. The 
bare bound state energy can be tuned by varying the 
magnetic field B 



E C (B) = 6n(B - B c ) 



(23) 



where 5/i = //atoms - /-tc is the magnetic moment difference 
and B c is the field where E C (B C ) = at threshold. The 
scattering length is real with 6 = and takes on the 
following resonant form 



a(B) 



abg — ab g 



B-Bq' 



where 



A = and B = B C + SB. 



(24) 



(25) 



Note that the interaction between the entrance and 
closed channels shifts the point of singularity of a(B) 
from B c to Bq . Such magnetically tunable Feshbach res- 
onances are characterized by four parameters, namely, 
the background scattering length ab g . the magnetic mo- 
ment difference 5fi, the resonance width A, and position 
B . 

Figure[B]shows an example of the scattering length and 
bound state energies for the 40 K 87 Rb molecule near the 
lowest energy spin channel of the separated atoms. The 
spin quantum numbers and hyperfinc splitting in their 
respective electronic ground states are 1, 2, and 6.835 
GHz for 87 Rb and 9/2, 7/2 and -1.286 GHz (inverted) 
for 40 K. There are 11 other closed spin channels in this 
system with E a > E\ that have the same total projec- 
tion quantum number as the lowest energy a = 1 spin 
channel. Because of their different magnetic moments the 
energy of a bound state of one of these closed channels 
can be tuned relative to the energy of the two separated 
atoms in the a — 1 s-wavc channel, as shown in the Fig- 
ure. Due to coupling terms in the Hamiltonian among 
the various channels, bound states that cross threshold 
couple to the entrance channel and give rise to resonance 
structure in its a(B). The resonance with Bo near 54.6 
mT (546 G) has been used to associate a cold 40 K atom 
and a cold 87 Rb atom to make a 40 K 87 Rb molecules in a 
near-threshold state with a small binding energy on the 
order of 1 MHz or less (l9j . 

It is extremely useful to introduce the properties of the 
long range van der waals potential and take advantage of 
the separation of short and long range physics discussed 
in the previous Section. Assuming that the interaction 
W(R) is confined to distances R <C a, the matrix element 
in Eq. (fT5|) defining T(E) can be factored as 



T(E) = C bg AE)- 2 t , 



(26) 
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FIG. 6: Molecular bound state energies (lower panel) and 
scattering length (upper panel) versus magnetic field B in mT 
(1 mT = 10 Gauss) for the lowest energy a — 1 s-wave spin 
channel of the 40 K 87 Rb fermionic molecule. The bound state 
energies are shown relative to the channel energy E\ of the two 
separated atoms taken to be zero. This a = 1 spin channel has 
respective 40 K and 87 Rb spin projection quantum numbers 
of —9/2 and +1, giving a total projection of —7/2. In this 
species there are 11 additional closed s-wave channels with 
E a > Ei and with the same projection of —7/2. The bound 
state quantum numbers are a(i), where i is the vibrational 
quantum number relative to the dissociation limit of closed 
channel a = 2, . . . , 12. Four bound states cross threshold in 
this range of B, giving rise to singularities in the scattering 
length. 



where f is a measure of resonance strength that depends 
only on the energy-independent short-range physics near 
E = 0, and is completely independent of the asymptotic 
boundary conditions. It thus can be used in character- 
izing the properties of both scattering and bound states 
when E =/= 0. 

The extrapolation of resonance properties away from 
E = depends on two additional parameters associated 
with the long range potential, jj, and Cq, which deter- 
mine a and E. Let us define a dimensionless resonance 
strength parameter 



aE 



(27) 



where r bg = a bg /a. Using the threshold van der Waals 
form of C bgi o(E)~ 1 given in the previous section, we can 
write 



1 



(28) 



The above-threshold scattering properties are found from 
the scattering phase shift rj{E) ~ r] bg (E)+ri res (E) 1 where 
r) res (E) is found from Eq. ([TTD once E c , T(E) and SE(E) 
are known. The first two are given by Eqs. (|2"3")l and (f2"tj)) . 



SE(E) = -t&n\ bg (E), 



(29) 



where tan X bg (E) is a function determined by the van 
der Waals potential, given a bg . It has the limiting form 
tanXbg(E) = 1 — r bg as E — > 0, and tanAb g (£') = for 
E ^> E [H, Q . Thus the position of the scattering length 
singularity is shifted by 



SB = B - B c = A 



i + (i-r bg y 



(30) 



from the crossing point B c of the "bare" bound state. 
Scattering phase shifts calculated from the van der Waals 
potential with the "quantum defect" forms in Eqs. (|26p 
and (|29[) are generally in excellent agreement with com- 
plete coupled channels methods for energy ranges on the 
order of E and even larger [111 ]. 

The properties of bound molecular states near thresh- 
old can also be calculated from the general coupled- 
channels quantum defect method using the properties of 
the long range potential. When the energy E b {B) = 
— h 2 kb(B) 2 /(2fi) of the threshold s-wave bound state is 
small, that is, \E b (B)\ < E or k b {B)a < 1, then the 
equation for E b {B) from the quantum defect method is 



(E C (B) - E b {B)) 



1 \ f 

- k b (B)a = - 

r bg -l J 2 



(31) 



If T = 0, we recover the uncoupled, or "bare," bound 
states of the system, whereas when T > 0, this equa- 
tion gives the coupled, or "dressed," bound states. The 
threshold bound state "disappears" into the continuum 
at B = B n , where a(B) has a singularity. The shift 
in Eq. |3"01 follows immediately upon solving for E c (Bq) 
where E b (B ) = 0. 

Threshold bound state properties are strongly affected 
by the magnitudes of s res and rb g . When the coupled 
bound state wavefunction is expanded as a mixture of 
closed and background channel components, |c) and \bg) 
respectively, an important property is the norm Z(B) 
of the closed channel component; the norm of the en- 
trance channel component is 1 — Z(B). The value of 
Z can be calculated from a knowledge of E b (B), since 
Z=\5 f i- 1 dE b /dB\ 0. 

There are two basic classes of resonances. One, for 
which s res 1, are called entrance channel dominated 
resonances. These have Z(B) <C 1 as B — Bo varies over 
a range that is a significant fraction of |A|. In addition, 
the bound state energy is given by Eq. (|11[) over a large 
part of this range. On the other hand, closed channel 
dominated resonances are those with s res <C 1. They 
have Z(B) large, on the order of unity, as \B — Bq\ varies 
over a large fraction of |A|, and only have a "univer- 
sal" bound state Eq. ([3]) over a quite small range <C |A| 
near Bq. Entrance channel dominated resonances have 
T(£', B) > E when < E < E, so that no sharp reso- 
nance feature persists above threshold, where a(B) < 
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FIG. 7: The lower panel shows an expanded view of Et(B) 
near B for the 40 K 87 Rb resonance with B = 54.693 mT 
(546.93 G) in Fig. [6] The solid line comes from a coupled 
channels calculation that includes all 12 channels with the 
same —7/2 projection quantum number. The dashed and 
dotted lines respectively show the universal energy of Eq. Q 
and the van der Waals corrected energy of Eq. . The up- 
per panel shows the closed channel norm Z(B). The width 
A = 0.310 mT (3.10 G), a bg = -191 a , and Sfi/h = 33.6 
MHz/mT (3.36 MHz/G). With a = 68.8 a and E/h = 13.9 
MHz, this is a marginal entrance channel dominated reso- 
nance with s re3 — 2.08. 

and the last bound state has disappeared. By contrast, 
closed channel dominated resonances with \rbg\ not too 
large will have T(E, B) < E when < E < E, so that 
a sharp resonant feature emerges just above threshold, 
continuing as a quasibound state with E > into the 
region where a(B) < 0. 

Figure [7] shows an expanded view of the 4(— 2) reso- 
nance of 40 K 87 Rb near 54.6 mT. The figure shows the 
character of the bound state as it merges into threshold 
at Bq. It tends to be a universal "halo" bound state 
over a range of \B — Bq\ that is less than about 1/3 of 
A. As \B — Bq\ increases, the bound state increasingly 
takes on the character of the closed channel 4(— 2) level 
as Z increases towards unity. Figure [H] shows an exam- 
ple of the very broad 6 Li resonance in the lowest energy 
a = 1 s-wave channel, which requires two 6 Li fermions 
in different spin states. This is a strongly entrance chan- 
nel dominated resonance, where Z < 1 over a range of 
\B — Bq\ nearly as large as A. The last bound state is 
a universal halo molecule over a range larger than 100 
G. The corrected Eq. (fTTj) is a good approximation over 
an even larger range. The scattering length graph shows 
that the size « a(B)/2 of the halo state is very large 
compared to a = 30 ao (see Table HJ over this range. 

Magnetically tunable scattering resonances have 
proven very useful in associating two cold atoms to make 
a molecule in the weakly bound states near threshold. 
This work is reviewed in detail in Ref. [l[ . The magneto- 




B (mT) 



FIG. 8: Molecular bound state energy (lower panel) and scat- 
tering length (upper panel) versus magnetic field B for the 
lowest energy a = 1 s-wave spin channel of the 6 Li2 molecule. 
This channel has one 6 Li atom in the lowest +1/2 projection 
state and the other in the lowest —1/2 projection state for 
a total projection of 0. There are 4 additional closed chan- 
nels with projection 0. In this range of B there is only one 
bound state that crosses threshold at Bo = 83.4 mT (834 G). 
The lower panel shows Eb(B) from a coupled channels cal- 
culation (solid circles) , the universal limit of Eq. ((3} (dashed 
line) and the corrected limit of Eq. (|ll|l (solid line). The 
width A = 30.0 mT (300 G), a bg = -1405 a , and Sfi/h = 28 
MHz/mT (2.8 MHz/G). This is a strongly entrance channel 
dominated resonance with s re3 = 59, and Z < 0.06 over the 
range of B shown. 

association process works by first preparing a gas with 
a mixture of both atomic species at B > Bq (assuming 
5(1 > 0), where there is no threshold bound state. By 
ramping the B field down in time so that B < Bq, col- 
liding pairs of atoms with E > can be converted to 
diatomic molecules in a bound state with energy E < 0. 
The conversion efficiency will depend on both the ramp 
rate and the phase space density of the initial gas. If 
the initial atom pair is held in a single cell of an opti- 
cal lattice instead of a gas, the conversion efficiency can 
approach 100 per cent. A simple Landau-Zener picture 
has been found to be quite accurate for such lattice cells, 
where the conversion probability of the atom pair in the 
trap ground state i = is 1 — e~ A , where 

2tt W 2 

A = — — ^ . (32) 
h E c 

Here i > represents the above-threshold levels of the 
atom pair confined by the trap, continuing the below 
threshold series of dimer levels with i < —1. For a 
three dimensional harmonic trap with frequency uj x = 
ujy = oj z = u>, the matrix element W C i = (C\W(R)\i) is 
well-approximated as W ci = y / T(E i ) /2ny / dE i /di, where 
dEi/di = 2hui and T(Ei) = 2fciah g <5/iA for the i = trap 
ground state of relative motion with k% = \fi^uj/h (sec 
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Eqs.[18l [2"0land [25]) . The trick used here in getting a ma- 
trix element W C i between two bound states from the ma- 
trix element (C\W(R)\E) involving an energy-normalized 
scattering state is to introduce the density of states as in 
Eq. pop . In a similar manner, the matrix element can 
be obtained between the bare closed channel state and 
the bound states i < of the entrance channel. Such 
matrix elements characterize avoided crossings like the 
one in Fig. [6]for E/h near —0.4 MHz and B near 43 mT. 
Finally, it should be noted that a Landau-Zener model 
can also be used for molecular dissociation by a fast mag- 
netic field ramp. An alternative phcnomcnological model 
has been developed to describe molecular association in 
cold gases, which are more complex than two atoms in a 
lattice cell El. 
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V. PHOTOASSOCIATION 

Cold atoms can also be coupled to molecular bound 
states through photoassociation (PA), as discussed in 
Chapters XX, YY, and ZZ. Fig. [9] gives a schematic de- 
scription of PA, a process by which the colliding atoms 
can be coupled to such bound state resonances through 
one- or two-photons. Reference Q reviews theoretical 
and experimental work on PA spectroscopy and molecule 
formation. Molecules made using the magnetically tun- 
able resonances described in the last Section are neces- 
sarily very weakly bound, with binding energies limited 
by the small range of magnetic tuning. Photoassocia- 
tion has the advantage that laser frequencies are widely 
tunable, so that a range of many bound states becomes 
accessible to optical methods, even the lowest v = vi- 
brational level of the ground state. In addition, the light 
can be turned off and on or varied in intensity for time- 
dependent manipulations. 

Photoassociation naturally lends itself to the reso- 
nant scattering treatment of a decaying resonance in 
Eq. (|2"2"j) . which applies to the one-color case with po- 
sition E c = E* — hvi, strength a b9 r (/) = T(E,I)/(2k), 
and shift SE(I). The latter two are linear in laser inten- 
sity / when / is low enough. PA is usually detected by 
the inelastic collisional loss of cold atoms it causes, due 
to the spontaneous decay of the excited state to make hot 
atoms or deeply bound molecules. In the limit E — ► 
the complex scattering length is 



a(vi,I) 



a bg 



J opt 



lcE 



El 



(7c/2) 2 



opt 



El 



(7c/2) 2 



(33) 
(34) 



where Eq = E* — hv\ + SE(I) is the detuning from res- 
onance, including the intensity-dependent shift, and the 
optical length is defined by L opt = ab g To(I)/j c . 

Photoassociation spectra, line shapes, and shifts have 
been widely studied for a variety of like and mixed alkali- 
metal species. At the higher temperatures often encoun- 



FIG. 9: Schematic representation of one- and two-color pho- 
toassociation (PA). The two colliding ground state atoms at 
energy E can absorb a laser photon of frequency u\ and be ex- 
cited to an excited molecular bound state at energy El. The 
bound state decays via spontaneous emission at rate 7 c /?i. If 
a second laser is present with frequency 1/2, the excited level 
can also be coupled to a ground state vibrational level v at 
energy E v , if h (i>2 — vi) = E — E v . The PA process depends 
on the ground state wavefunction at the Condon point Rc of 
the transition, where hvi equals the difference between the 
excited and ground state potentials. 



tered in magneto-optical traps, contributions to PA spec- 
tra from higher partial waves, e.g., p- or d- waves, have 
been observed in a number of cases. The theory can be 
readily extended to higher partial waves. By introducing 
an energy-dependent complex scattering length the the- 
ory for s-waves can be extended to finite E away from 
threshold and to account for effects due to reduced di- 
mensional confinement in optical lattices [2lj |. 

The optical length formulation of resonance strength 
is very useful for a decaying resonance. It also applies 
to decaying magnetically tunable resonances, if Tq from 
Eq. ([25]) is used to define a resonance length ab g S/jA/j c 



equivalent to L op t [22|. The scattering length has its 
maximum variation of ab g ± L opt when the laser is tuned 
to Eq = ±7 c /2, and losses are maximum at Eq = 
where b = L opt . When detuning is small, on the or- 
der of 7c, significant changes to the scattering length 
on the order of a are thus normally accompanied by 
large loss rates (see Eq. [14]) . Losses can be avoided 
by going to large detuning, since when {^ c /Eq) <C 1, 
b = (L opt /2)('j c / 'Eq) 2 , whereas the change in a only varies 
as a — abg = —L opt {^ c /EQ). To make the change a — a,b g 
large enough while requiring {^ c /Eq) <C 1 means that 
L op t has to be very large compared to a. 

The magnitude of L opt depends on the matrix ele- 
ment (C\h,fli(R)\E) where hfli(R) represents the optical 
coupling between the ground and excited state. Using 
Eqs. ([TH]) and ([2D)) and the above definition of L opt , and 
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factoring out the relatively constant TiVt\ value, 



L op t = n 



l^il 2 F(E) 



lc k 

The Franck-Condon overlap factor is 

poo 

F{E) = / ip*(R)tfj (R,E)dR 
Jo 



(35) 

(36) 
(37) 



where Dq is the derivative of the difference between 
the excited and ground state potentials evaluated at the 
Condon point Rc, and dE*/dv is the excited state vi- 
brational spacing. Equation (|37|) is known as the re- 
flection approximation, generally an excellent approx- 
imation where F(E) is proportional to the square of 
the ground state wavefunction at Rc, the Condon point 
where the molecular potential difference matches hv\ (see 
Fig. [5]). Thus, F(E) can be evaluated using expressions 
like Eqs. J5| or j9f for Rc 3> a or Rc <C a respec- 
tively. The reflection approximation is quite good over 
a wide range of E and for higher partial waves than the 
s-wave. By selecting a range of excited levels v by chang- 
ing laser frequency u±, thus changing Rc, the shape and 
nodal structure of the ground state wavefunction can be 
mapped out over a range of R. 

The optical length has several important properties ev- 
ident from Eq. (|35[) . First, since both fli and 7 C are 
proportional to the same squared transition dipole mo- 
ment, L opt does not depend on whether the transition 
is strong or weak, but can be large for both kinds of 
transitions. Second, L opt oc |f2i| 2 so it can be increased 
by increasing laser intensity. Third, since F{E) cx k as 
E — ► for entrance channel s-waves, L opt oc F(E)/k is 
independent of E or k at low energy. However, it does 
depend strongly on the molecular structure through the 
Franck-Condon factor. In practice, using strong transi- 
tions with large decay rates such as those in alkali-metal 
species leads to the requirement to use excited molecu- 
lar levels far from threshold with large binding energies. 
This is necessary to achieve large detuning from atomic 
and molecular resonance. This requirement means such 
levels have small F(E) factors, due to the very large value 
of Dc in Eq. {37}. On the other hand, weak transitions 
with small decay rates, such as those associated with the 
So -^ 3 Pi intercombination line transition for alkaline 
earth species such as Sr, can lead to quite large values of 
L opt . This is because large detuning in 7 C units can be 
achieved for levels that are still quite close to the excited 
state threshold. Such levels typically have large Franck- 
Condon factors. In fact, PA transitions near the weak in- 
tercombination line of Sr have been observed to have L opt 



several orders of magnitude larger than was observed for 
strongly allowed molecular transitions involving Rb [23[ . 
Thus, there are good prospects for some degree of optical 
resonant control of collisions in ultracold gases of species 
like Ca, Sr, or Yb. 

Two-color PA is also possible when a second laser with 
frequency v-i is added, as shown in Fig. [5J When the 
the frequency difference is chosen so that h{v2 — v\) = 
E—En, the ground ii molecular level is in resonance with 
collisions at energy E. By keeping v\ fixed and tuning 
v-2, two-color PA spectroscopy can be used to probe the 
level energies. This is how the data on binding energies 
of the Y\>2 molecule was obtained [l5[ so as to be able 
to construct Fig. In this case, 12 different levels from 
different isotopic species were measured, among which 
were levels with i = — 1 and —2 and I = and 2. Two 
color spectroscopy has also been carried out for several 
alkali-metal homonuclear species. 

Two-color processes are also an excellent way to assem- 
ble two cold atoms into a translationally cold molecule. 
Early work along these lines was done using the sponta- 
neous decay of the excited level to populate a wide range 
of levels in the ground state. The disadvantage of spon- 
taneous decay is that it is not selective. However, by 
using a laser with a precise frequency, a specific level can 
be chosen as the target level. One early experiment did 
this to associate two 87 Rb atoms in a Bosc-Einstein con- 
densate to make a molecular level at a specific energy of 
ft(-636) MHz [13. 

It is highly desirable to be able to make translationally 
cold molecules in their vibrational ground state v = 0. 
This is especially true of polar molecules, which have 
large dipole moments in v — 0. On the other hand, 
threshold levels have negligible dipole moments, since 
there is no charge transfer because of the large average 
atomic separation ps a. A promising technique is to use 
magnctoassociation using a tunable Fcshbach resonance 
to associate the atoms into a threshold molecular level, 
then use a 2-color Raman process to move the popula- 
tion in that state to a much more deeply bound level. 
Although molecules in a gas are subject to fast destruc- 
tive collisions with cold atoms or other molecules in the 
gas (see Eq. IT?)) , the molecules can be protected against 
such collisions by forming them in individual optical lat- 
tice trapping cells. Then the 2-color Raman process could 
be used to produce much more deeply bound molecules 
that are stable against destructive collisions. This has 
been done successfully with 87 Rb2 25j molecules. In the 
future, such methods are likely to produce v = polar 
molecules, with which a range of interesting physics can 
be explored [H, [27} . 
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